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Method for converting digital image pixel values 



(57) A method is provided for converting values rep- 
resentative of pixel intensities from an input dynamic 
range to an output dynamic range. The input dynamic 
range is determined by acquisition circuitry (20). while 
the output range is determined by display or output de- 
vices or circuitry (36). A histogram of the input values is 
generated, and a lower limit value for the useful portion 
of the input range is identified from the histogram. An 



upper limit of the useful input range is then identified 
based on a log transformed histogram computed from 
the original histogram, and on a threshold value located 
based on the lower limit value. Conversion of the^ pixel 
values Is then performed by reference to the lower and 
upper limit values. The technique facilitates compensa- 
tion for variations in utilization of the input dynamic 
range, and appropriately scales the input range utilized 
to the available output range. 
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Description 

[0001] The present invention relates to digital imaging 
and, more particularly, to the conversion of a range of 
input pixel values to a different range of output pixel val- 
ues representative of an image, such as an X-ray image. 
[0002] Discrete imaging devices, such as digital X-ray 
imaging systems, employ a detector which divides re- 
gions of an image into individual picture elements, or 
pixels. The array or matrix of pixels defines, when 
viewed as an overall image, features of interest, such 
as internal anatomy of a subject positioned adjacent to 
the detector. To facilitate interpretation by physicians 
and technicians, the individual intensities of the pixels 
typically define the features of the image by Imitating 
contrasts and textures obtainable through conventional 
film-based X-ray or imaging systems. 
[0003] To convert the detected pixel intensities to dig- 
itized values suitable for display, the pixel intensity val- 
ues are processed after acquisition by the detector. In 
a first stage, the detected pixel intensities are digitized 
in values which vary over a predetermined dynamic 
range of the detector and acquisition circuitry, such as 
12 to 14 bits. In X-ray systems, for example, these dig- 
itized values are representative of the quantity of X-rays 
received by each pixel during data acquisition. Subse- 
quently, the pixel intensity values are scaled to map the 
values onto the dynamic range of a display device. As 
part of this scaling, it is common to perform logarithmic 
transformation of the image pixel values to obtain a re- 
sulting image which mimics conventional film-rendered 
images. In addition, the scaling process maps the dy- 
namic range of the detector and acquisition circuitry on- 
to the dynamic range of the display. The latter range may 
be substantially different from that of the upstream cir- 
cuitry, such as on the order of 8 to 10 bits. 
[0004] While the logarithmic transformation of the dig- 
itized pixel values is useful in rendering an image which 
is understandable by attending physicians and techni- 
cians, performing the transformation prior to the dynam- 
ic range scaling can be problematic. For example, his- 
tograms are often employed to analyze pixel intensity 
values. However, processing of histograms generated 
based on the transformed values can result in difficulties 
in identifying high and low limits of relevant portions of 
the detected data, rendering the dynamic range scaling 
difficult. The use of logarithmically transformed data pri- 
or to dynamic range scaling can also result in loss of 
accuracy for individual pixels in the image matrix. 
[0005] There is a need, therefore, for an improved 
method for processing discrete image data that facili- 
tates use of as much of the dynamic ranges of acquisi- 
tion circuitry and display circuitry as possible. In partic- 
ular, there is a need for an improved method for convert- 
ing digital pixel values defining a discrete pixel image 
from a first dynamic range to a second dynamic range 
in a computationally efficient manner. 
[0006] In an exemplary emt>odimenf of the invention, 



digitized values for a plurality of pixels of an image, such 
as an X-ray image, are converted from a first dynamic 
range to a second dynamic range for eventual display 
of the overall image. The input values, distributed over 

5 the first dynamic range, are used to form a histogram 
representative of the number of pixels having predeter- 
mined digital intensity values. A lower limit of the rele- 
vant portion of the input data dynamic range is identified 
from the histogram. A transformed histogram is devel- 

10 oped based upon the Input value histogram. The trans- 
formed value histogram is used to determine a threshold 
value for the log-transformed image. Based upon this 
threshold value, an upper limit of the relevant portion of 
the input data range is identified from the transformed 

IS value histogram. With the lower and upper limits of the 
input value dynamic range thus Identified, the input pixel 
values are nnapped onto output values over the second 
dynamic range. 

[0007] The technique offers a computationally effi- 
20 cient approach to mapping of input and output dynamic 
ranges in image processing, particularly in digital X-ray 
systems. By identifying the low and high limits of the in- 
put dynamic range, the useful range of the input values 
is identified and the mapping is performed over this use- 
2S ful range. The technique may be employed for individual 
images, thereby accommodating variations in intensity 
values for individual patients and individual Images. The 
resulting images provide consistent appearance, facili- 
tating comparison and interpretation of imaged features. 
30 [0008] The invention will now be described in greater 
detail, by way of example, with reference to the draw- 
ings, in which :- 



Fig. 1 is a diagrammatical representation of an X- 
35 ray imaging system including circuitry for converting 
input values for individual pixels of an image into 
output values for display; 

Fig. 2 is a flow chart illustrating steps and exemplary 
control logic for converting or scaling digitized pixel 

40 values over a first dynamic range to values over a 
second dynamic range for display; 
Fig. 3 is an exemplary histogram of pixels having 
predetermined intensity values for an image detect- 
ed via the system of Fig. I in accordance with the 

45 logic of Fig. 2; and 

Fig. 4 is an exemplary histogram of transformed pix- 
el values generated based on the histogram of Fig. 
3, in accordance with the logic summarized in Fig. 2. 

so [0009] Ref emng now to Fig. 1 , an imaging system 10 
is illustrated in the fonm of an X-ray system. It should be 
noted, however, that the technique disclosed herein is 
not limited to use in such X-ray systems, but may find 
application in other types of digitized image processing 

55 modalities. Imaging system 10 includes a X-ray source 
12 configured to emit a stream of X-rays 14. When an 
image is to be captured, a subject, such as a human 
body 16 is positioned between the X-ray source 12 and 
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an image acquisition unit 16. Image acquisition unit 18 
includes a detector 20 and an intensity encoding circuit 
22. A portion of the X-rays emitted by the source 1 2 pass 
through the body and strike the detector 20. 
[0010] Detector 20 is configured to divide the surface 
receiving the X-rays into a matrix of discrete picture el- 
ements or pixels 24. Detector 20 outputs a stream of 
data signals representative of the intensity of radiation 
received by each of the discrete pixel locations. Intensity 
encoding circuit 22 receives the data signals from de- 
tector 20 and encodes the intensity levels in a data 
stream in which the intensity level or value of each pixel 
is assigned a digitized value over a predetermined dy- 
namic range. For example, intensity encoding circuit 22 
may digitize the intensity values over a dynamic range 
of 12 to 14 bits per pixel, providing 2^^-2^^ different dig- 
itized intensity values potentially assignable to each pix- 
el. 

[0011] Signal receiving unit 18 then transmits these 
encoded values to an Image processing unit 26 for fur- 
ther refinement. Image processing unit 26 Includes cir- 
cuitry configured to process the digitized values and to 
render output data which can be displayed in a form 
comprehensible to an atteriding physician or technician. 
In the diagrammatical illustratbn of Fig. 1, Image 
processing unit 26 includes a decoding circuit 28, a cor- 
rection circuit 30, a value conversion circuit 32. and a 
signal processing circuit 34. Decoding circuit 28 re- 
ceives the digitized values from intensity encoding cir- 
cuit 22 and partially decodes the information, such as 
by placing the pixel values in a desired order or se- 
quence, as will be appreciated by those skilled in the 
art. Circuit 28 then transmits the decoded values to cor- 
rection circuit 32 which performs corrections on the val- 
ues, such as to produce digital intensity values which 
are linearly proportional to the quantity of X-rays re- 
ceived in each pixel location of detector 20 during the 
acquisition process. As will be appreciated by those 
skilled in the art, correction circuit 30 may perform var- 
ious corrections on the values, such as to normalize the 
values to produce, optimum, consistent presentations 
of images acquired by the system, thereby compensat- 
ing for such factors as the size of the patient, the com- 
position of the tissues or structures of the body, the X- 
ray dosage, and so forth. 

[0012] The corrected digital intensity values are then 
transmitted to value conversion circuit 32. Value conver- 
sion circuit 32 converts the digital intensity values from 
the range output by intensity encoding circuit 22 to a 
range appropriate for display, as summarized below. 
These converted values are then further processed by 
signal processing circuit 34, which may perform such 
functions as enhancing structural details, smoothing 
texture, and so forth. Signal processing circuit 34 then 
outputs values for the digitized image to a display/output 
device 36, typically including a monitor or printer capa- 
ble of producing a discrete pixel Image 38 representa- 
tive of the features of interest within bcSdy 16. 



[001 3] Value conversion circuit 32 accommodates dif- 
ferences in the dynamic ranges of values output by sig- 
nal receiving unit 18 and values needed by display/out- 
put device 36. For example, while pixel intensity values 
5 output by unit 18 may have a dynamic range of 12-14 
bits per pixel, the input range of display/output device 
36 may be smaller, on the order of 8-10 bits per pixel. 
To permit both acquisition and display circuitry to utilize 
the greatest portion of their respective dynamic ranges, 
10 then, value conversion circuit 32 analyzes the useful 
portion of the dynamic range of the Incoming digitized 
values, and, based upon this useful portion, converts the 
incoming digitized values to the dynamic range required 
by display/output device 36. 
IS [0014] It should be noted that certain of the functional 
circuitry illustrated in Fig. 1, including correction circuit 
30 and value correction circuit 32, may be embodied in 
appropriate code within a general purpose or applica- 
tion-specific microprocessor or computer. Moreover, the 
microprocessor or computer may serve to carry out oth- 
er functions and signal processing operations beyond 
those described herein. 

[0015] Fig. 2 represents steps in exemplary control 
logic implemented by value conversion circuit 32 in de- 
termining the useful portion of the incoming digitized im- 
age pixel values and in converting or scaling these val- 
ues to the dynamic range of downstream circuitry. As 
Illustrated In Fig. 2, once value conversion circuit 32 re- 
ceives the digitized signals representative of pixel inten- 
sities, a histogram h(x) is generated to determlne.pixel 
populations having individual Intensity values over the 
full input dynamic range. A schematic illustration of a 
histogram of this type Is illustrated in Fig. 3. 
[0016] As shown in Fig. 3. the histogram 52 may be 
represented graphically by plotting intensity values 54 
across a horizontal axis, and population counts 56 for 
each of the intensity values along a vertical axis. The 
resulting graphical representation forms bars 58 which 
define a distribution curve 60 in a stepwise manner. 
Each bar 58 has a width W such that the entire dynamic 
range is defined over the horizontal axis between a dig- 
itized value of zero bits to a maximum value of 2' - 1 bits, 
where "i" represents the number of bits in the input dy- 
namic range. It should also be noted that the sum of the 
population counts for all of the intensity values equals 
the total population of the pixels considered In the im- 
age. 

[0017] While certain images may include pixels hav- 
ing Intensities over most of the useful Input range (0 to 
2^ - 1), value conversion circuit 32 analyzes the input 
values to determine the useful portion of the dynamic 
range over which the values lie. Thus, where variations 
In the acquisition process result in utilization of less than 
all of the input dynamic range, the converted values will 
nevertheless provide as much information as is availa- 
ble from the input values. 

[0018] In a first step in this process, as indicated by 
step 42 in Fig. 2, circuit 32 determines the lower limit L 
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of this useful range. This step is graphically represented 
in Fig. 3, In Fig. 3. reference numeral 62 indicates an 
exemplary lower limit L established by reference to his- 
togram 52. The lower limit value is selected by locating 
an intensity value along the horizontal axis below which s 
a desired percentage of the global pixel population falls. 
That is, beginning with an intensity value of zero, the 
pixel counts 56 of bars 58 for each intensity value are 
accumulated until the desired pixel population percent- 
age is obtained. Once the percentage Is reached, the 
lower limit is set at the next higher intensity value. In 
general, lower limit values corresponding to 0.5 to 1 per- 
cent of the global pixel population have been found to 
provide satisfactory results. 

[0019] Referring again to Fig. 2, once the lower limit 
L has been identified, a log-transformed histogram h(f 
(x)), defined as g(x), is calculated as indicated at step 
44. An exemplary log-transformed histogram 64 is illus- 
trated in Fig. 4. The log-transformed histogram may be 
generated from the equation: 

h(f(x)) = xh(x): 



value, such as 0.8T, and the foregoing procedure is per- 
formed using this reduced threshold to identify a value 
72 further along the histogram 64. 
[0022] As indicated by step 48 in Fig. 2, based upon 
the value 72, an upper limit value H for the useful portion 
of the input dynamic range is next determined. The val- 
ue of H is set to a multiple of the value 72 determined 
as summarized above. It should be noted that this mul- 
tiple need not be an integer multiple, and will, in general, 
differ for various types of modalities and subjects. For 
example, in digital X-ray imaging systems, the multiple 
of the value 72 used to set the upper limit H will differ in 
accordance with a need to view the skin/air interface, 
lung interstitium, or other details. By way of example, in 
chest radiographs, a multiple of 1 .5 has provided good 
results over a range of images. 
[0023] With the lower and upper limits of the useful 
portion of the input dynamic range thus determined, val- 
ue conversion circuit 32 next scales the pixel values 
from the input range to the output range as indicated at 
step 50 in Fig. 2. In the presently preferred technique 
this conversion is performed in accordance with the fol- 
lowing relationship: 



15 



where the values of x and h(x) are based upon the val- 
ues of the pixel intensity histogram 52 of Fig. 3. The log- 
transformed histogram may be smoothed using a low- 
pass filter technique, such as by local averaging of the 
population counts for each intensity value with the im- 
mediately preceding and succeeding population counts. 
Such filtering tends to smooth peaks and valleys which 
may occur in the histogram profile. 
[0020] Once the log-transformed histogram g(x) has 
been established, a threshold value T is determined, as 
indicated at step 46 in Fig. 2, Referring again to Fig. 4, 
the threshold value T is located from the log-trans- 
formed histogram 64 by identifying a value of g(x), as 
indicated by reference numeral 66, at which the lower 
limit value T, indicated by numeral 62, intersects the log- 
transformed histogram, as indicated at Intersection 
point 68 in Fig. 4. The histogram 64 is then scanned from 
the value x = L until a second intersection point 70 is 
located at which the histogram value g(x) falls below the 
threshold value T This value will form the basis for es- 
tablishing an upper limit for the useful portion of the input 
dynamic range. 

[0021] Before establishing the upper limit value, how- 
ever, a verification comparison is preferably performed 
to determine whether a satisfactory portion of the dy- 
namic range has been selected through use of the 
threshold value T. In particular a comparison of the iden- 
tified value 72 and a multiple of the lower limit value has 
been found to provide a satisfactory test of the width of 
the dynamic range. For example, if the value identified 
at 72 is determined to be less than 2L, the algorithm may 
be considered to have found a spurious feature of the 
histogram, requiring use of a different threshold value 
T In such cases, the threshold value is reset to a lower 



where the value Output is the scaled output value for 
each pixel, Input \s the corresponding input value for the 
30 pixel, b is the number of bits in the output dynamic range, 
e is small number (less than or equal to 1 .0, and typically 
0.5) used to assure that the output value does not go 
below zero, and L and H are the lower and upper limit 
values, respectively. The brackets represent the opera- 
35 tion of (determining the largest integer value less than or 
equal to the real-valued quantity enclosed therein. Note 
that the transformation described by this equation re- 
verses high and low pixel values, such that pixels re- 
ceiving low input values will have high output values (re- 
40 suiting in inversion of light and dark pixels in the result- 
ing image). This transformation is appropriate, for ex- 
ample, in images such as chest radiographs. However, 
as will be appreciated by those skilled in the art, other 
transformations which do not reverse high and low val- 
45 ues may be appropriate in other applications. 

[0024] It has been found that the foregoing method 
provides satisfactory conversion of input values to out- 
put values in a computationally efficient manner. More- 
over, the procedure is readily adaptable to various im- 
50 aging modalities, as well at to a wide variety of input and 
output circuitry having substantially different dynamk; 
ranges. Finally, the method automatically accommo- 
dates variations in the bandwidth of the dynamic range 
actually utilized, even in particular images produced by 
55 a single imaging system or on a particular subject, pro- 
viding more consistent and comparable images, and 
thereby facilitating image interpretation and analysis. 
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Claims 

1. A method for converting digitized image pixel val- 
ues from a first range to a second range, the method 
comprising the steps of.: 

(a) determining a lower limit value of a relevant 
portion of the first range; 

(b) generating an intensity histogram repre- 
sentative of pixel populations having specified 
intensities, and transforming the histogram to 
generate a log-transformed histogram; 

(c) identifying a threshold value for an upper 
limit of log-transformed values from the log- 
transformed histogram; 

(d) identifying a population of pixels having log- 
transformed values having a desired relation- 
ship to the threshold value; 

(e) determining an upper limit value of the rel- 
evant portion of the first range based upon the 
Identified population; and 

(f) converting the pixel values to converted val- 
ues over the second range based upon the low- 
er and upper limit values. 

2. A method for determining a useful range of image 
pixel intensity values in a digital pixel imaging sys- 
tem, the method comprising the steps of: 

(a) determining a bwer limit value of the useful 
range by identifying an intensity value below 
which values of a desired fractional portion of 
the pixels fall; 

(b) generating an intensity histogram and a log- 
transformed histogram based upon the intensi- 
ty histogram; 

(c) identifying a threshold value for an upper 
limit value based upon the log-transformed his- 
togram; 

(d) Identifying a population of pixels having log- 
transformed values having a desired relation- 
ship to the threshold value; and 

(e) determining an upper limit value of the use- 
ful range based upon the identified population. 

3. The method of claim 1 or 2, wherein the population 
identified in step (d) is a population having log- 
transformed values lower than the threshold value. 

4. The method of claim 1 , 2 or 3, wherein the lower 
limit value is determined based upon populations of 
pixels identified in the intensity histogram. 

5. The method of any preceding claim, wherein the up- 
per limit value is determined based upon the log- 
transformed histogram and the population Identified 
in step (d). 



6. The method of any preceding claim, wherein the 
threshold value corresponds to a desired log-trans- 
formed value. 

s 7, The method of any preceding claim, wherein the 
population Identified in step (d) is identified by com- 
parison of the threshold value to log-transformed 
values in the log-transformed histogram. 

10 a The method of any preceding claim, including the 
further steps of comparing the population identified 
In step (d) to a desired population size, and resetting 
the threshold value based upon the comparison. 

IS 9. The method of claim 8, wherein the desired popu- 
lation size is determined based upon a population 
of pixels having a predetermined relationship with 
the lower limit value. 

20 10. A method for converting digitized image pixel val- 
ues from an input dynamic range to an output dy- 
namic range, the method comprising the steps of: 

(a) generating a first histogram of the pixel val- 
25 ues; 

(b) identifying a lower limit value from the first 
histogram; 

(c) generating a second histogram of log-trans^ 
formed values of the pixel values by transfor- 

30 nnation of the first histogram; 

(d) identifying a population of pixels having de- 
sired log-transformed values from the second 
histogram; 

(e) identifying an upper limit value from the first 
3S histogram based upon the identified popula- 
tion; and 

(f) converting the pixel values from the input 
range to converted values over the output 
range based upon the lower and upper limit val- 

40 ues. 

11. The method of claim 10, wherein the population is 
identified in step (d) by reference to a desired 
threshold log-transformed value. 

46 

12. The method of claim 10 or 11, comprising the further 
step of comparing the identified output range to a 
desired output range, and modifying the desired 
threshold log-transfonned value based upon the 

so comparison. 

13. The method of claim 12, wherein the desired output 
range is determined based upon the lower limit val- 
ue. 

55 . 

14. The method of any one of claims 10 to 13, wherein 
the lower limit value Is identified based upon a value 
below which values of a desired fractional portion 
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(54) Method for converting digital image pixel values 



(57) A method is provided for converting values rep- 
resentative of pixel intensities from an input dynamic 
range to an output dynamic range. The input dynamic 
range is determined by acquisition circuitry (20), while 
the output range is determined by display or output de- 
vices or circuitry (36). A histogram of the input values is 
generated, and a lower limit value for the useful portion 
of the input range is identified from the histogram. An 



upper limit of the useful input range is then identified 
based on a log transformed histogram computed from 
the original histogram, and on a threshold value located 
based on the lower limit value. Conversion of the pixel 
values is then performed by reference to the lower and 
upper limit values. The technique facilitates compensa- 
tion for variations in utilization of the input dynamic 
range, and appropriately scales the input range utilized 
to the available output range. 
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